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Abstract 



The basic features of the slow relaxation phenomenology arising in phase 
ordering processes are obtained analytically in the large iV model through the 
exact separation of the order parameter into the sum of thermal and conden- 
sation components. The aging contribution in the response function Xag(t, t w ) 
is found to obey a pattern of behavior, under variation of dimensionality, qual- 
itatively similar to the one observed in Ising systems. There exists a critical 
dimensionality (d = 4) above which Xag(t,t w ) is proportional to the defect 
density po(t), while for d < 4 it vanishes more slowly than pD(t) and at 
d = 2 does not vanish. As in the Ising case, this behavior can be understood 
in terms of the dependence on dimensionality of the interplay between the 
defect density and the effective response associated to a single defect. 
PACS: 05.70.Ln, 64.60. Cn, 05.40.-a, 05.50. +q 

I. INTRODUCTION 

The phase ordering processes |lj following the quench of non disorder systems (e.g. fer- 
romagnets) below the critical point provide a simplified framework for the study of slow 
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relaxation phenomena. In particular, aging and the off equilibrium deviation from the fluc- 
tuation dissipation theorem (FDT) have been studied numerically |2|-|J and through analyt- 
ical treatment of solvable models The basic structure of the complex phenomenology 
arising in these processes originates from the wide separation in the time scales of fast and 
slow variables. Referring to the more intuitive case of a domain forming system, in the late 
stage of phase ordering the order parameter can be assumed to be the sum of two statistically 
independent components f| 

(f)(x,t) = ip{x,t) + a(x,t) (1) 

the first describing equilibrium thermal fluctuation within domains and the second off- 
equilibrium fluctuations due to interface motion. From ([J) follows rather straightforwardly 
the split of the autocorrelation function 

G(t,t w ) = G st (t-t w )+G a JpJ (2) 

where G s t(t — t w ) is the stationary time translation invariant (TTI) contribution due to i/j 
and G ag {t/t w ) the aging contribution due to the off-equilibrium a degrees of freedom. A 
similar structure shows up in the linear response at the time t to an external random field 
switched on at the earlier time t w 

X(t,t w ) = Xst(t-t w ) +Xag(t,t w ). (3) 

Here, the stationary contribution Xst{t — t w ) satisfies the equilibrium FDT with respect to 
G st (t — t w ), while Xag{t,t w ) is the off-equilibrium extra response due to the presence of the 
interfaces. This is considered [@,^]|| to be proportional to the interface density 

Xag(t,t w ) ~ pj(t) (4) 

where pi(t) ~ L^it) and Lit) ~ t 1 ^ is the typical domain size with z = 2 for non conserved 
order parameter |IJ, as we shall assume in the following. A more precise formulation of this 
behavior is through the scaling relation 
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Xa g (t,t w )=tZ a X {£) (5) 

with the exponent a = 1/2. The implication is that Xag(t,t w ) is negligible in the asymptotic 
regime t w — > oo. 

Motivated by analytical results |7||§ for the one dimensional Ising model which do not 
fit in the above scheme, recently we have undertaken a detailed study of the behavior of 
the response function under variation of space dimensionality for a system with a scalar 
order parameter j|. On the basis of numerical simulations for discrete Ising spins and 
approximated analytical results for continuous spins, we have arrived to a picture for the 
behavior of Xag(t,t w ) which modifies considerably the one presented above. This is best 
understood by introducing the notion of the effective response due to a single interface and 
defined by 

Xag(t, t w ) = pi(t)Xeff(t, t w ). (6) 

Then, if (||) were to hold, Xeff(t,t w ) ought to be a constant. Instead, we have found that 
for an Ising system this is the case only for d > 3, while for d < 3 there is the power law 
growth 

Xeff(t,t w )~(t-t w ) a (7) 

with numerical values for the exponent compatible with a = (3 — d)/4. At d = 3 the 
power law is replaced by logarithmic growth, promoting d = 3 to the role of a critical 
dimensionality. This result is interesting for two reasons. The first concerns the mechanism 
of the response and the role of the off-equilibrium degrees of freedom. The power law (0) 
reveals that for d < 3 the aging component of the response does not originate trivially 
just from the polarization of interfacial spins, but a more complex phenomenon producing 
large scale optimization of the position of domains with respect to the external field is at 
work. The second regards the overall behavior of Xag(t,t w ). Putting together (||) and (|7|) 
the exponent a in (H) acquires a dependence on dimensionality 
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showing that lim^^oo Xag(t, t w ) does not vanish as d — > 1. This indicates that d = 1 plays the 
role of a lower critical dimensionality, where the off equilibrium response becomes persistent. 
An interesting consequence of this phenomenon is that the connection between static and 
dynamic properties [|Kj], which holds for d > 1, is invalidated at d = 1. It should be added 
that in order to have a phase ordering process in d = 1 thermal fluctuations within domains 
must be suppressed [|J. 

As stated previously, the picture summarized above has been established on the basis 
of a combination of exact results for the d = 1 Ising model, numerical results for Ising 
systems with d > 1 and approximate analytical results for continuous spin systems. It is, 
then, interesting to test how general is the picture. As a step in this direction, we have 



considered the large N model where exact analytical calculations can be carried out [11 



The slow relaxation properties of the large N or equivalent mean field models, arising in the 
quench at or below Tc have been analyzed before jp|,|r2|]. What we do here, however, goes 
beyond previous results since we manage to reproduce exactly and analytically the scenario 
outlined for Ising systems. We show that in the large N model one can make explicitly the 
separation ([]]) of the order parameter into the sum of two independent components which are 
responsible of the stationary and aging contributions in (fj). Then, we carry out analytically 
the corresponding separation ([3D of the response function. After introducing the notion of 
defect density for the large N model, we derive by analogy with (||) the effective response 
per defect and we find a behavior, as dimensionality is varied, qualitatively similar to the 
one established for scalar systems. 

The paper is organized as follows. In Section [TT] the large N model is defined and the 
main static properties are reviewed. In Section [TTI| the solution of the equation of motion 
is presented and the analytical form of the autocorrelation function is obtained in quenches 
at and below the critical point. In Section [TV] the splitting of the order parameter into 
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independent components satisfying the requirements above described is carried out. In 
Section [V] the behavior of the integrated response function is considered in relation to the 
off equilibrium deviation from the FDT. Finally, concluding remarks are made in Section [VT|. 



II. MODEL AND STATIC PROPERTIES 

We consider the purely relaxational dynamics of a system with a non conserved order 
parameter governed by the Langevin equation 

^P = -3^ + r/(f,t) (9) 
dt S(f>{x,t) 

where <fi = (0i, </>jv) is an iV-component vector, ff(x,t) is a gaussian white noise with 
expectations 

(f)(x,t)) = 

(10) 

( Va (x , t) V/3 (x', f)) = 2T5 a ^6(x - x')5(t - f ) 

and T is the temperature of the thermal bath. In the dynamical process of interest the 
system is initially prepared in the infinite temperature equilibrium state with expectations 

(f(x)) = 

m (11) 

(<j) a {x)(j)j3{x')) = A5 a ^8(x - x 1 ) 

and at the time t = is quenched to a lower final temperature Tp. The hamiltonian is of 
the Ginzburg-Landau form 



n[(j)}= / d d x 

J \f 



I( V 0) 2 + -0 2 + ^-(^) 2 



(12) 



where r < 0, g > and V is the volume of the system. In the large N limit the equation of 
motion for the Fourier transform of the order parameter <p{k) = J v d d x(p(x) exp(ik ■ x) takes 
the linear form 

= -[k 2 + I(t)]$(k,t)+ff(k,t) (13) 

where 



(ff(k,t)) = 

WV 77 (14) 
(Va(k,t) V p(k',t')) = 2T F 8 a , /3 V5 %+ p Q 5(t - f) 

and the function of time 

/(t)=r + |(0 2 (f,t)> (15) 

must be determined self-consistently, with the average on the right hand side taken both over 
thermal noise and initial condition. If the volume V is kept finite the system equilibrates in 
a finite time t eq and the order parameter probability distribution reaches the Gibbs state 

P eq [}$)\ = I e -^£^-^>>« (16) 

z 

where £ is the correlation length defined by the equilibrium value of I(t) through 

r 2 = r + |(0 2 (x)) e9 (17) 

with (-) eg standing for the average taken with flllf) . 

In order to analyze the properties of P eq [<p(k)] it is necessary to extract from fllTD the 
dependence of on T and V. Evaluating the average, the above equation yields 

k 



The solution of this equation is well known [T3|] and here we summarize the main features. 



Separating the k = term under the sum, for very large volume we may rewrite 

C 2 = r + gT F B(C 2 ) + g^z- 2 (19) 

where 

_ fc 2 

9x ,. 1 1 f d d k e'^ 1 , nn . 

B ^-^v^WT^ = J(^?WT^ (20) 

regularizing the integral by introducing the high momentum cutoff A. The function B{x) is 
a non negative monotonically decreasing function with the maximum value at x = 
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™=J$n?- = w- i ih A " ■ (21) 

By graphical analysis one can easily show that ([T^) admits a finite solution for all Tp. 
However, there exists the critical value of the temperature Tc defined by 

r + gT c B(0) = (22) 

such that for T F > T c the solution is independent of the volume, while for T F < T c it 
depends on the volume. Using 

B{x) = (^rixl-^T (l - i, ^ (23) 

where r(l — |, -jfe) is the incomplete gamma function, for < Tf f ^ c <S 1 one finds (Appendix 
I) £ ~ ( Tf tJ P )~ U wnere v — 1/2 for d > 4 and z/ = l/(d - 2) for d < 4, with logarithmic 
corrections for d — 4. At Tc one has £ ~ U A with A = 1/4 for d > 4 and X — l/d for c? < 4, 
again with logarithmic corrections for d = 4. Finally, below Tc one finds £ 2 = where 
M 2 = M 2 ( S ^ k ) and M 2 = -r/#. 

Let us now see what are the implications for the equilibrium state. As ([16]) shows the 
individual Fourier components are independent random variables, gaussianly distributed 
with zero average. The variance is given by 

^($(k)-$(-k)) eq = VC eq {k) (24) 

where 

CeM = p^z, (25) 

is the equilibrium structure factor. For Tp > Tc, all k modes behave in the same way, with 
the variance growing linearly with the volume. For Tp < Tc, instead, £ -2 is negligible with 
respect to k 2 except at k — 0, yielding 

2g(l - 6 {n ) + aV 2X 5r n , for T F = T c 
C eq {k) = { kK k ' 0) M (26) 

W~ko)+ M2V ko ,forT F <T c 



where a is a constant (Appendix I). This produces a volume dependence in the variance 
of the k = mode growing faster than linear. Therefore, for Tp < Tc the k = mode 
behaves differently from all the other modes with k ^ 0. For Tp < Tc the probability 
distribution ( |T6| ) takes the form 

P^m-^-S^e-^^^^ . (27) 

Zj 

Therefore, crossing Tc there is a transition from the usual disordered high temperature phase 
to a low temperature phase characterized by a macroscopic variance in the distribution of 
the k = mode. The distinction between this phase and the mixture of pure states, obtained 
below Tc when N is kept finite, has been discussed elsewhere [O . 



We shall refer to this transition as condensation of fluctuations in the k = mode. In 
order to gain a better insight it is convenient to go back to real space, splitting the order 
parameter into the sum of the two independent components 

0(f) =a + $(x) (28) 



with 



and 



Then (|27 ) takes the form 



where 



a = ±$(k = 0) (29) 



*)4^i e ' fe ( 30 ) 



P e9 [0(f)] = P{a)P[^{x)\ (31) 



P^) = TK^FKWH^ ( 32 ) 



(2vrM 2 ) iV / 2 



shows the formation of the condensate below Tc with the macroscopic variance ^(c? 2 ) e <? 
M 2 while 



p[^(f)] = I e -w/v^(V^ (33) 

Zj 

describes thermal fluctuations about the condensate. Correspondingly, the correlation func- 
tion G eq {x — x') = (1/N)(<p(x) ■ 4>{x')) eq splits into the sum of two pieces 

G eq (x -x ) = G T (x - x') + M 2 (34) 

where 

G T {x-x') = ^{x)-^x')) eq (35) 

is the correlation of thermal fluctuations and M 2 comes from fluctuations of the condensate. 
Since Gt(x-x') at large distances decays like \x— x'\ 2 ~ d , from (|34]) follows lim^^^^ G eq {x— 
x') = M 2 showing the violation of the clustering property of the correlation function due to 
the breaking of ergodicity in the low temperature phase. For future reference, notice that 
from flTT| ) follows that for x = x! and Tp < Tq 

G eq (0) = M 2 (36) 

which in turn implies 

G T (0) = M 2 - M 2 . (37) 

As stated above, with a finite V equilibrium is reached for t ~ t eq and t eq ~ Hence, 
in a quench to Tp < Tc one has t eq ~ V implying that if the V — > oo limit is taken 
at the beginning of the quench, equilibrium is not reached for any finite time. Since in 
the following we are interested in the relaxation regime before equilibrium is reached, we 
shall take the thermodynamic limit from the outset. We are interested to see whether it 
is possible to carry out the decomposition ([I]) of the order parameter in such a way that 
a(x, t) eventually evolves into the equilibrium condensate a and ip(x, t) into the equilibrium 
thermal fluctuations. 
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III. DYNAMICS 



Due to rotational symmetry, from now on we shall drop vectors and refer to the generic 
component of the order parameter. The formal solution of (|i"3]) is given by 



</>(k, t) = R(k, t, 0)M%) + / dt'R(k, t, t')r](k, f) (38) 

JO 

where 4>o{k) = <f)(k,t = 0) and according to (pi]) in the V — » oo limit 

The response function is given by 

R(k,t,t') = ^e- k2 ^ (40) 

with Y(t) = exp[Q(t)], Q(t) = Jq dsl(s) and y(0) = I. The actual solution is obtained once 
the function Y(t) is determined. In order to do this, notice that from the definition of Y(t) 
follows 

dY 2 (t) r i 

^P = 2[r + g(<p\x,t))}Y\t). (41) 

Writing (<p 2 (x,t)) in terms of the structure factor 

d d k 



<^fc*)> = J j^y d c ( k ^ e ~ p (42) 

and using ( |3~8|) to evaluate C(k,t) 

C{k, t) = R 2 {k, t, 0)A + ITp J dt'R 2 {k, t, t') (43) 
from fl4~l|) we obtain the integro-differential equation 

^ = 2rV(«) + 2 S A/ (« + ^L) + 4 9 T F /J *-/ + ^) K»(0 (44) 



where 



f{x) = I = (Bkx)-*. (45) 
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Solving (f44|) by Laplace transform [[T^,|12|], the leading behavior of Y(t) for large time is 
given by 

A a e^' H ; for T F > T c 

Ar ■ for Tp = T c (46) 
A b t~i ; for T F < T c 

where u = for d > 4 and u = (d — 4)/2 for d < 4. The values of the constants A a , Ab and 
yl c are listed in Appendix II. 

This completes the solution of the model. Once the response function is known, we may 
go back to (|38D and take various averages. Using (|TJ) and ( |39D we have 



(0(M)> = o 



(47) 



and for the two time structure factor (<f)(k,t)<f)(k' ,t')) = C(k, t, t')(27r) d 5(k + k') with t > t' 



C{k, t, t') = R{k, t, 0)R(k, t', 0)A + 2T F / dt"R(k, t, t")R{k, t', t"). 



(48) 

The corresponding real space correlation function G(x — x',t,t') = (4>(x, t)4>(x', t')) is given 
by 

G(x - x', t, t') = J dx"R(x - x", t, 0)R(x' - x", t', 0) A 

+ 2Tp J* dt" J dx"R(x - x", t, t")R(x' - x" , t', t") (49) 

where R(x,t,t') is the inverse Fourier transform of R(k,t,t'). In the following we will be 
primarily concerned with the autocorrelation function G(t,t') = G(x — x! = 0,t,t'). Using 
the definitions fl4H|) and (^5|) this is given by 



G(M') 



1 



Y(t)Y(t') 



t + t' 



Y 2 (f" 



2A 2 



(50) 



The behavior of this quantity for different final temperatures and for different time regimes 



has been studied in the literature 15,121. Here we summarize the results. 



For Tp > Tc from ( Pf ) follows that for t' > t eq = 2£ 2 the autocorrelation function is 



TTI 
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G(r) = G eq (Q)e 



(51) 



where G eq (0) is the equilibrium fluctuation above Tc given by (|17D , i.e. G eq (0) = Mq + ^£~ L> 
and t = t — t'. 

For Tp < T c the equilibration time diverges and there are two time regimes of interest: 
i) short time separation: t' — > oo , £ — > 
u) large time separation: t' — ► oo , £ — > oo. 



Taking i' large and using (fK3), for Tp = Tc in these limits we find 



G(t,t') 



M 2 (A 2 r + l) 1 "! , for ^ -> 
At /1 -iT(^) , for £ -> oo 



with 

F(ar) = < 
and for Tp < T c 

G(t,t' 



(47r) d /2(d-2) 
4 

(87r) d / 2 (d-2) 



(¥) 



!- rf / 2 /'s+lN 1 -' 1 / 2 ' 



for 2 < d < 4 
for d > 4 



M 2 + (M 2 - M 2 )(A 2 r + I) 1 -* , for £ -»• 



M 2 



i/f 



.(*+? ) 2 

In the latter case we may also write 



, for — > oo. 



G(t,t') = G st (r) + G ag ^jj 



where 



G st (r) = (M 2 -M 2 )(A 2 r + l) 1 ^ 



and 



G ag (x) = M 2 



Ax 



1 + x) 



(52) 



(53) 



(54) 



(55) 



(56) 



(57) 



This is illustrated in Fig.[l| which shows the convergence toward the form fl5"5| ) of the exact 
autocorrelation function G(t, t') obtained by solving numerically the coupled set of equa- 



tions (H) and (|g). 
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As explained in the Introduction, this behavior is suggestive of the existence of two 
variables responsible respectively of the stationary and of the aging behaviors. We want 
now to show that in the large N limit these two variables can be explicitly constructed. 

IV. SPLITTING OF THE FIELD 

The task stated at the end of the previous Section requires the splitting of the order 
parameter field into the sum of two independent contributions 

<t>(x,t) = ip(x,t) + a(x,t) (58) 

with zero averages (ip(x,t)) = cr(x,t) = and autocorrelation functions such that 

(i>(x,t)il)(x,t')) = G st (T) (59) 

a(x,t)<j(x,t') = G ag (t'/t) . (60) 

In order to stress the statistical independence of the two component fields, we have used the 
angular brackets for averages over ip and the overbar for averages over a. 

For the construction of fields with these properties, let us go back to Fourier space. Using 
the multiplicative property of the response function 

R(k, t, t')R(k, t', t ) = R(k, t, t ) (61) 

with t > t' > to it is easy to show that the formal solution (|38|) of the equation of motion can 
be rewritten as the sum of two statistically independent components <j)(k, t) = ip(k } t)+a(k, t) 
with 

a(k,t) = R(k,t,t )(/)(k,t ) (62) 

and 

i/)(k,t)= /' dt'R(k,t,t')r](k,t') (63) 
13 
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since for < t < t, (p(k,t ) and 7](k,t) are independent by causality. In other words, the 



order parameter at the time t is split into the sum of a component a(k,t) driven by the 
fluctuations of the order parameter at the earlier time to and a component ip(k,t) driven by 
the thermal history between to and t. Let us remark that to can be chosen arbitrarily between 
the initial time of the quench (t = 0) and the observation time t. With the particular choice 



ip component describes fluctuations of thermal origin while the a component, as it will be 
clear below, if t is chosen sufficiently large, describes the local condensation of the order 
parameter. 

According to definitions (g2j) and (H), from (0) and © follows a(k,t) = (ip{k,t)) = 0. 
The two time structure factor splits into the sum 



t = 0, the component a(k,t) is driven by the fluctuations in the initial condition (|39|). The 



C(k, t, t') = C a (k, t, t') + C^k, t, t') 



(64) 



with a(k, t)a(k', t') = C a (k, t, t') (2n) d 5(k + k') and 



C a (k, t, t') = R(k } t, t Q )R(k, t', t )C(k, t ) 



(65) 



where C(k, to) is the equal time structure factor at the time to. Similarly (ip(k, t)ip(k', t')) 
C 1 p(k,t,t')(2iT) d 5(k + k') with 




(66) 



Going to real space and setting x = x', we have 



G(t,t') = G a (t,t')+G4t,t') 



(67) 



with 




(68) 



and 




(69) 
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Assuming that t and t' are sufficiently larger than t so that C(k,t ) under the integral 
in (Pf) can be replaced by its value at k = and using (EH) we may write 



G a {t,t') 



Y 2 {t ) Jt + t' I 



Y(t)Y(t 



a/ 



2A 2 



(70) 



Gf(t, t) 



2A 2 



(71) 



where C = C{k = 0,t ). 

Let us now evaluate these results for short and large time separations. In the first case 
the dominant contributions are given by 



G a (t,t') = M 2 — (M 2 - M 2 )(2AV^ 



(72) 



and 



GV(M') = (M 2 - M 2 ) [(2A%) 1 ^ + (A 2 r + l) 1 ^ 



(73) 

where the unknown constant Cq entering (ffOl) has been eliminated imposing that the equi- 
librium sum rule G a (0) + G^(0) = Mq be satisfied. The sum of the above contributions is 
independent of t , as it should, and coincides with (|54l) . Similarly, in the large time regime 
we find the t dependent results 



G a (t,t') = \M 2 - (Ml - M 2 ){2A 2 t y-^ 



4tf 



(t + t>) 



(74) 



and 



GV(M') = (M 2 - M 2 )(2A 2 t ) 1 ^ 



At? 



(75) 



(t + t>)\ 

again with the sum independent of to and giving back (|54]) . Defining the microscopic time 
t* = A -2 , we see that taking t 3> t* from ([F2|) and (|73|) we get the short time behavior 



G a (t,?) = M 2 

G4t,t') = (M 2 -M 2 )(A 2 T + l) 1 -l 



(76) 
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(77) 



Comparing with ( f)E| ) and (|57| ) we can make the identifications 



G^t') = G st (r) 



(78) 



and 



G a (t, t') — G, 



ag 



( 



) 



(79) 



Therefore, the fields a(x, t) and ^(af, t) defined by (|62D and (|63D , with the choice of i ^ £*, 
provide an explicit realization of the decomposition ([58]) satisfying the requirements ( |59"|) 
and (|60|). The physical meaning of the two components can be readily understood comparing 
the equal time values of (|76|), which yield respectively G a {t = t') = M 2 and G^p(t = t') = 
Mq — M 2 , with the equilibrium results ([36]) and (^). It is then clear that the field a(x, t) is 
associated to local condensation of the order parameter with fluctuations of size M 2 , while 
the field ip(x,t) describes thermal fluctuations. In the case of a system with a scalar order 
parameter, a would be associated to the average value of the order parameter characteristic 
of a domain and ip to thermal fluctuations within domains ||. This makes also clear the 
origin of the requirement to 3> t*. In order to make a separation of variables with the above 
physical meaning, it is necessary to wait a time t large enough for all microscopic transients 
to have occurred leaving well formed local equilibrium. 



In the previous Section we have produced the explicit separation of the order parameter 
into the condensation component and thermal fluctuations in the quench at Tp < Tc- It 
is now interesting to see in what relation these components are with the linear response 
function. 



V. RESPONSE FUNCTION 
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If an external field is switched on at the time t w > to, the splitting fl58|) modifies into 
4>h{x, t) = ip(x, t) + (i h (x, t) where to linear order 

a h (x, t) = a(x, t) + C dt' f dx*R(x - x', t, t')h(x') (80) 

with a(x,t),ij;(x,t) and R(x,t,t') unperturbed quantities. Here, we are interested in the 
response to a quenched, gaussianly distributed random field with expectations 

E h [h(x)} = 

L 1 ' l (81) 
E h [h a (x)hp(x')} = hl5 aP 5{x - x' ). 

Computing the staggered magnetization from (BO) and averaging over the field we obtain 

N ^2y J y dxE h [a h (x, t) ■ h{x)\ = x(t, t w ) (82) 

where x(t^w) — It w ^(^^') * s the integrated response function and R(t,tf) = R(x — x! = 
0,t,t'). From ((p and fl4|) for T F < T c 

_d d 

, /x f d d k ,-> .. fc 2 . , d ft'\ 4 / , 1\"2 

R(t, = J j^y R ( k > *> * ) e ~^ = ^y~ z ( 7 ) v ~ av • (83) 

Let us then write the integrated response function as the sum 

X(t,t w ) = Xst(t - t w ) + Xag(t,t w ) (84) 

where the stationary component Xst(t — t w ) is defined by requiring that the equilibrium FDT 
be satisfied with respect to the stationary component ( |78| ) of the autocorrelation function, 
namely 

T FX st{t-t w ) = G^)-G i) {t-t w ). (85) 

It is straightforward to check that this is verified by 

Xst{t - t w ) = (4tt)-2 jf dt' (t-t'+ — ' 

= ^« - M 2 ) {l - [A 2 (t - *„) + l] 1 ^} . (86) 

The aging component then remains defined by the difference 
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d 



Xag(t, t w ) = x(t, t w ) - Xst(t - t w ) 

= {A-n)~n w 2 x~ l / dy(l-y 



x 



A 2 t 



(87) 



where x = t w /t. This shows that for d > 2 and any fixed value of x lim tTO _ i . 0O Xas(£> — 
implying 



lim = - 



(88) 



Hence, in the limit t w — > oo the equilibrium FDT (|85|) is satisfied by the whole response 
function and the plot of x(t,t w ) vs G^(t,t w ) is linear. Instead, if x(t,t w ) is plotted against 
the full autocorrelation function (0), from (|55|) follows 



lim T F x{t,t v 

t w — >OQ 



G(t, t) - G(t, t w ) , for M 2 < G(t, t w ) < M 2 



M 2 - M 2 



(89) 



, for G(t,t w ) < M 2 

yielding the behavior of Fig.|2| which is characteristic of the phase ordering process [fj,[| . 

For d = 2 the power of t w in front of the integral disappairs from (5^) and the leading 
contribution for t ^> A~ 2 is given by 



(90) 



showing that the aging contribution to the response does not vanish as t w — > oo. Further- 
more, since T c = for d = 2, the phase-ordering process requires T F = and from 
follows 



G(t,t')=G a (t,t') = 2M 2 
Eliminating x between ( p0|) and ([31]) we get 

1 



x 



1 + x 



(91) 



Xa S (G) 



2tt 



log 



1 + 



G 



1 - 



G 2 



(92) 



which gives a parametric plot (Fig. |3|) qualitatively similar to what one finds in the Ising 
model for d — 1. 



As it was recalled in the Introduction, with a scalar order parameter the behavior of 
Xag{t, t w ) under variations of dimensionality becomes trasparent by introducing the effective 
response Xeff(t,t w ) associated to a single interface. The mechanism regulating the behavior 
of Xag(t,t w ) then is explained through the balance between the rate of loss of interfaces as 
coarsening proceeds, and the rate of growth of the single interface response given by (^). 
The dimensionality dependence of the exponent a in that case is the outcome @] of the 
competition between the external field and the curvature of interfaces in the drive of interface 
motion. 

In the large N model there are no localized defects and none of the above concepts has 
direct physical meaning. Nonetheless, the notion of defect density can be extended to the 
large N case by looking at the behavior of I(t). Consider the quench to Tp = where there 
are no thermal fluctuations. From (|46"D, for large time, I(t) ~ —t^ 1 namely 

M 2 -±(fi(x,t)}~t-\ (93) 

If the system was in the ordered state, with the order parameter aligned everywhere, the 



left hand side ought to vanish. Therefore, the positive difference (|93|) may be attributed to 
"defects" with density 

p D (t) ~ L~ 2 (t) ~r l (94) 

which is what one obtains in general with a vector order parameter |IJ . We may then define 
the effective response function per defect by the analogue of (|6]) Xa g (t,t w ) = p D (t)x e ff{t, t w ) 
which yields 

i ~i 1 

Xeff(t,t w )=t 2 -% I dy j. (95) 

The behavior of this quantity can be computed analytically for short and for large time 
separation obtaining in both cases, apart from a change in the prefactor, a power law 
behavior as in (|^) 

X eff(t,t w )~(t-t w ) a (96) 
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with 



) , for d > 4 
a = { (97) 

2 - \ , for d < 4 
and 

Xeff(t,t w ) ~ \og [A 2 (t-t w )} (98) 

for c? = 4. The overall exact behavior of Xeff(t,t w ) is depicted in Fig.£|, obtained by plot- 
ting ( P^D for different dimensionalities. 

Therefore, the qualitative picture is the same obtained in the scalar case [|]]. There exist 
upper critical dimensionalities, d\j = 3 for N — 1 and dj/ = 4 for N = oo (presumably for 
all iV > 1), above which Xeff saturates to a constant value within microscopic times. At djj 
Xeff grows logarithmically, while below dy there is power law growth. In addition there are 
lower critical dimensionalities, di, — 1 for N — 1 and di = 2 for iV = oo (or N > 1), where 
the exponent a reaches exactly the value such that Xeff(t,t w ) ~ Pj 1 ^) in the scalar case 
and Xeff{t,t w ) ~ Pd (t) in the vectorial case. Namely, at dt the growth of Xeff makes up 
exactly for the loss of interfaces or defects, producing the behavior illustrated in Fig. [| 



VI. CONCLUDING REMARKS 

In this paper we have shown that the basic features of the slow relaxation phenomenology 
arising in phase ordering processes: separation of order parameter, autocorrelation function 
and linear response function into fast and slow components, can be obtained analytically 
and exactly in the large N model. The behavior of Xag{t, t w ) is of particular interest since it 
displays the same qualitative pattern of behavior under variation of dimensionality observed 
in the Ising case. This is a strong indication that this might be a generic feature of the 
slow relaxation in phase ordering processes. In this respect it might be worth to undertake a 
numerical study of Xag{t, t w ) in systems with vector order parameter and finite N at different 
dimensionalities. 
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A comment should be made about the connection between static and dynamic properties. 
One of the most interesting recent developments in the study of the off equilibrium deviation 
from FDT has been the derivation [|Hj of a link between the response function and the 
structure of the equilibrium state. Assuming that lim^^oo t w ) = x(G(t,t w )), i.e. that 
in the large time regime x{t, t w ) depends on time only through the autocorrelation function, 
this connection takes the form 

/ G=q 

where P(q) is the probability that in the equilibrium state the overlap -^y f v dx<f)(x) ■ <f>'(x) 
between two different configurations [0(x)] and [0'(^)] fakes the value q. For Ising sys- 
tems holds for d > 1. In fact, for d > 1 the t w — > oo limit of x(t,t w ) is found [@,f|f§ to 
have the form (|89|) which, applying (P^), yields 

P(q) = 8(q - M 2 ) (100) 

in agreement with an equilibrium state formed by the mixture of pure states. For d = 1, 
since Xag(t,t w ) does not disappear as t w — > oo, the connection (|99|) is no more valid ||. 
Then, in the large iV model we should expect (|99D to fail at most for d = 2, since that is the 
case where Xag(t,t w ) does not asymptotically vanish. However, it is not so straightforward 
that (|99|) should hold for d > 2. In fact, as observed above, the behavior of x(t,t w ) in the 
limit t w — >• oo for d > 2 is indistinguishable from what one obtains in the Ising case. Namely, 
one finds the form (p9|) of x(G) which yields ( |100j ) for P(q) and this is what one expects 
when there is a mixture of ordered pure states. The problem is, as explained in Section 2, 
that the structure of the low temperature state in the large iV model is quite different from 
a mixture of pure states. This puzzling feature of the large iV model will be investigated in 
a separate paper. 
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VII. APPENDIX I 



Solving (|T9|) in the large volume limit one finds: 



for < ^^b. < 1 



and 



(£A) 



-2 



A! 

2 



A/ 2 (87T) d / 2 

2T F r(l-d/2) 



e j f 



T F -Tg 1_ 



T C T(2-d/2) 



I 7b 



2 

" 2-d 



2 

d-2 



T F 2 

T C d-4 



, for d < 2 
, for d = 2 
, for 2 < d < 4 
, for d > 4 



(^A)- 2 log((eA)- 2 ) 



7>-T c 



for d — 4. 



For Tp = T c 



and 



T C \d-4 r y V 

Y(^)T^ l Mlk 2 - d 



for d > 4 



d , for d < 4 



£[21og(£A)]- 



A 2 T, 



for d = 4. 



VIII. APPENDIX II 



The pref actors in fl4"E|) are given by 



A n 



A 



8tt \ 

r(i-d/2) V?W 



M 2 



d/2' 



2 

' 2-d 



27r^-M 2 e- 4 ^ 

^-Tc)] |E| (i + ^) 

^ fj_ , AM£\ j 
2T C \2g ^ 2T C J log(2£" 2 ) 



Ml 

T C 



2T F |r(l-d/2)| 

(87T) d / 2 



d-2 1 

d-2 



V2g 2T C / V29 A 2 T c d-4 



-1 



, for d < 2 

, for d = 2 
, for 2 < d < 4 
, for d = 4 
, for d > 4 
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and 



A c = < 



(8n) d / 2 

A 



, for d < 2 
, for d = 2 



sm 



A = 



T F + 9 AM 2 



(8tt)<<, [M| (l - 



2 ' 



(106) 



(107) 
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FIGURES 




FIG. 1. Comparison of the exact autocorrelation function G(t,t w ) (broken line) with the sum 
G s t(t — t w ) + G ag (t/t w ) (continous line) for t w = 1, 10, 10 2 , 10 3 increasing from left to right. For 
t w = 10 3 the two curves are indistinguishable. Parameters of the quench are d = 3,Tp = Tc/2 
and A = 1. 
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FIG. 4. Plot of Xeff(t,t w ) for t w = 10 8 and A = 1. The dashed lines are power laws with the 
corresponding exponent a. 
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